Spin, charge and orbital ordering in Lao.sSri 5Mn04 
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We have analyzed the experimental evidence of charge and orbital ordering in Lao.5Sri.5Mn04 using 
first principles band structure calculations. Our results suggest the presence of two types of Mn 
sites in the system. One of the Mn sites behaves like an Mn^^ ion, favoring a Jahn- Teller distortion 
of the surrounding oxygen atoms, while the distortion around the other is not a simple breathing 
mode kind. Band structure effects are found to dominate the experimental spectrum for orbital and 
charge ordering, providing an alternate explanation for the experimentally observed results. 
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In recent times there has been a resurgence of inter- 
est in the manganites as a result of the wide range of 
physical properties that they exhibit. In some of the 
doped compounds, there is a real space ordering of the 
charge carriers in certain orbitals, resulting in orbital or- 
dering (00) and sometimes charge ordering (CO). As 
the magnetic and transport properties are closely cor- 
related with the orbital and charge degrees of freedom, 
there have been intense efforts - both experimental and 
theoretical - to clarify the microscopic mechanisms giv- 
ing rise to such exotic phenomena . While techniques 
such as electron microscopy and neutron scattering have 
been traditionally used to detect CO as well as mag- 
netic ordering the direct observation of 00 has been 
difficult. Murakami et al. performed pioneering experi- 
ments on LaMnOa |^ and Lao.5Sri.5Mn04 j|] and were 
successfully able to detect 00. The technique they used 
exploited the anisotropy of the x-ray scattering tensor at 
an absorption edge, thereby making the intensity of the 
scattered x-rays sensitive to the valence electron distri- 
bution. 

Lao.5Sri.5Mn04 is a layered manganite with an aver- 
age valence of 3.5 at each Mn site. In the high temper- 
ature phase, above the CO temperature, all Mn-0 bond 
lengths within the ab plane are equal [|6|. Neutron scat- 
tering studies by Sternlieb et al. ^ indicated a real space 
ordering of two distinct Mn sites below 217 K, with the 
two sites identified as Mn'^^ and Mn^+ sites. The inten- 
sity variation observed at the two sites was interpreted 
as arising from a 1% breathing mode- like movement of 
the oxygen atoms in the ab plane towards one set of Mn 
atoms, identified as the higher valent Mn''^+ sites. Later 
experiments by Murakami et al. used the anisotropy 
of the x-ray scattering tensor to show that the Mn sites 
form a pattern of alternating Mn^+-Mn^+ sites in the ab 
plane, with both CO and 00 taking place at the same 
temperature, 217 K. They claimed that the charge mod- 
ulation between the two Mn sites is not small, but corre- 
sponds to as much as an integral charge fluctuation be- 
tween neighboring Mn sites. However, if this is the case, 
strong correlations are necessary to localize charge car- 
riers. Surprisingly, however, certain other spectroscopic 



properties of this system can be easily explained by a 
single particle approach. 

The contrasting experimental evidence where one set 
of spectroscopic data could be easily explained by an ef- 
fective single particle approach, while another set im- 
plied the importance of strong correlations in the system 
prompted us to perform a detailed analysis of the elec- 
tronic structure using the density-functional pseudopo- 
tential framework. Apart from providing a microscopic 
mechanism for the observed orbital and charge-ordering 
in this system, our results have important implications 
on recent anomalous x-ray scattering experiments on 
Pri_,Ca,Mn03 @, LaSrzMnzOy |] and NaVzOs |o|. 

The details of the computational method are described 
elsewhere and we only mention minimum additional 
aspects. The cutoff energy for the plane wave expansion 
of the eigen state is 20 Ry, being smaller than the value 
of 30.25 Ry in our usual calculations because of the large 
unit cell in the present calculation. The k-point mesh 
used is 4x4x2. The generalized gradient approximation 
(GGA) was used for the exchange functional |12|. The 
alloying effects of (Lao.sSri.s) at the A-site of the lattice 
are treated within the virtual crystal approximation [ p^ . 

The Mn spins order below 110 K forming one- 
dimensional ferromagnetic zigzag chains in the ab plane 
coupled antiferromagnetically to each other, known as 
the CE-type magnetic ordering. As the exchange split- 
ting is large (about 3 eV), the amplitude for electrons 
to hop from a zig-zag chain of parallel spins (indicated 
by thick solid lines in Fig. la) to the neighboring chain is 
small. Hence the electronic structure is governed by these 
one-dimensional chains. Considering a single chain, it 
was shown ||l^ that the anisotropic hopping between the 
Mn d orbitals could explain the 00 at the two Mn sites. 
We have performed ab-initio calculations considering the 
complete 3-dimensional Immm structure [||. Similar to 
the findings in ref. we also found the kind of 00 that 
is observed experimentally. Choosing the x and y axes as 
shown in Fig. 1, we found that while the c?32.2_^2 orbital 
was preferentially occupied on atom 1, the d^y2_^-2 orbital 
was preferentially occupied on atom 3. These results can 
be understood within the framework of a simple nearest 
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neighbour tight-binding model involving the eg orbitals: 
d3a;2_j.2 and dj^2_^2 on atom 1, d3y2_j,2 and (1^2 on 
atom 3 as well as d^z^-r^ and dx2-y2 on atoms 2 and 4 
as used in Ref. [Q. While the d^x2_^2 orbital at Mn(l) 
and d-iy2_j.2 orbital on Mn(3) hybridize with both Cg or- 
bitals on atoms 2 and 4, the dy2_z2 on Mn(l) and dz2_^2 
on Mn(3) do not. The eigenvalue spectrum of such a 
simple model, consists of two bonding bands, energeti- 
cally separated from four non-bonding bands, with the 
antibonding bands at higher energies. As the average 
valence of Mn is 3.5, there are two electrons in the Cg 
orbitals of the 4 Mn atoms comprising the chain. Hence 
the two bonding bands which have dominantly d^x'^_j.2 
character on atom 1 and d-^y2_^2 character on atom 3 
are occupied. Thus, the one-dimensional chains which 
are a consequence of the magnetic structure drive the 
00 within these calculations. However, the difficulty 
in approaching the problem this way is that the Immm 
structure exists only above the CO/00 temperature and 
there are extensive evidence in the literature suggesting a 
symmetry lowering below the ordering transition. How- 
ever, the details of the crystal structure of this compound 
is shrouded in controversy. It was suggested jsj that the 
structure, within the charge ordered state, is Cmmm with 
lattice constants equal to a, \/2 a, c where a and c 
are the cell dimensions of the basic Immm cell. More re- 
cent work by Larochelle however, indicate that the 
structure is Ammm or one of its two subgroups Am2m or 
A222, with the lattice parameters being -\/2a, 2\/2 a, and 
c. None of the structural investigations have been able to 
provide the atom positions, particularly for the oxygen 
sites, owing to the large unit cell and the existing data 
quality [|l6| . This work however suggests certain possible 
distortions of the oxygens, though not in agreement with 
ref. In view of such uncertainities, we have optimised 
the structure of this system using ab-initio band struc- 
ture calculations. The structure that we find (see Fig. 
1) is P2i/m with lattice parameters \/2a, 2-\/2a and c. 
The lattice parameters as well as the space group that 
we have obtained are the same as in ref. JT6t, if the vir- 
tually negligible {'--^ 0.015 A) displacements of the Mn'^^ 
atoms shown in Fig 1 are ignored. 

On optimizing the internal coordinates, the system ex- 
hibited different distortions of the MnOg octahedra as- 
sociated with different Mn atoms. The direction of dis- 
placement of the oxygen atoms in the ab-plane is indi- 
cated by the arrows shown in Fig. lb, though the oxygen 
atoms have been left out of the figure for added clar- 
ity. The Mn site, Mn(3), showing d3y2_j.2 00 lowered 
its energy by an elongation of the Mn-0 bonds in the y- 
direction. This kind of distortion can be understood eas- 
ily within the framework of crystal- field effects. A tetrag- 
onal distortion of a MnOg octahedron resulting in an 
elongation of the Mn-0 bondlengths in the y-direction, 
lowers the bare energy of the d-^y2_^2 orbital. Thus, the 
sites Mn(l) and Mn(3) behave like Mn'^^ species, sustain- 
ing a Jahn- Teller (JT) like distortion of the surrounding 



oxygens in two mutually perpendicular directions, giving 
rise to the orbital ordering. As a consequence, the oxygen 
atoms surrounding the sites labelled Mn(2) and Mn(4) 
sustain distortions with the oxygen atoms along the x- 
axis being displaced in one direction and the ones along 
the y-axis being displaced in another direction, as shown 
in the figure, in contrast to the suggestion in ref. |^ where 
all four oxygen atoms surrounding the Mn*+ site move 
closer to the Mn atom. We found the kind of distortion 
suggested in that work to have a higher energy. The sites 
Mn(2) and Mn(4) have been identified as Mn^"*" sites in 
the literature and we retain that nomenclature, although 
within our calculations the charge difference between the 
so-called Mn^^^ and Mn^'^ species is negligible. 

While our calculations suggest different distortions of 
the surrounding oxygens about the Mn^+ and Mn'^^ sites, 
the actual magnitude of the distortion is fairly small . 
The underestimation of the magnitude of the JT distor- 
tion by these first-principles approaches is well-known, 
though these methods do get the nature of distortions 
correctly. An example of this is LaMnOs where the na- 
ture of distortion was correctly predicted, while the mag- 
nitude of the theoretical JT distortion was found to be 
half of the experimental value We simulated the 

neutron diffraction pattern using the optimized coordi- 
nates as well as the coordinates assuming a Jahn- Teller 
distortion equal to what is found in LaMnOa. The simu- 
lated diffraction patterns were nearly identical, indicating 
the difficulty of determining the oxygen positions with 
any precision from such measurements. Therefore the 
enhanced distortion of the oxygen octahedra are not in- 
consistent with the structural data available so far. We 
in fact hope that the present theoretical considerations 
will help in the future structural investigations. It is also 
significant that with optimized JT distortion, the correct 
ground state magnetic order cannot be reproduced for 
LaMnOa [|ll| . The situation seems to be the same in the 
present case, where the CE type antiferromagnetic state 
is more stable than the ferromagnetic state only with the 
enhanced JT distortion, consistent with earlier observa- 
tions [0. Hence in the subsequent analysis, we have 
extrapolated the magnitude of the JT distortion to the 
value observed in LaMnOa. 

The experiments of Murakami et al. [Q analyzed the 
energy dependence of the 00 superlattice reflection 
across the Mn K absorption edge. The strong energy de- 
pendence was interpreted as arising from a splitting of ~ 
5 eV between the Ap^ and 4py PDOS at the sites Mn(l) 
and Mn(3) which show 00. While the Ap^ states were 
suggested to be 5 eV lower than the Apy on one Mn atom, 
the order was reversed on the other Mn atom. Theoret- 
ically there are two contrasting interpretations possible 
for such a splitting. Ishihara and Maekawa |jl^ argued 
in the context LaMnOa that the splitting in the Mn Ap 
states induced by the intraatomic p~d Coulomb interac- 
tion produces such a strong tensor character of the scat- 
tering form factor. However, it was pointed out ||2^ that 
as the 4p states are extended, the suggested p-d Coulomb 
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interaction strength is unphysically large and the spht- 
ting in the 4p states is actually caused by the JT distor- 
tion, though Coulomb interactions, particularly within 
the Mn 3d manifold, may have important consequences 
for other properties in such systems. These suggestions 
are consistent with an earlier analysis of transition metal 
K-edge XAS in Fe and Co oxides ^ . 

In Fig. 2 we plot the Mn Ap^ and 4py partial densities 
of states (DOS) in the energy range 10-20 eV above the 
Fermi level, projected onto the Mn site (Mn(l)) showing 
dsx^-r^ ordering. A splitting of ~ 3 eV is clearly visible, 
with the Apx states located at lower energies. The order 
is reversed at the other Mn site (Mn(3)) which shows 
d^y2_j.2 ordering. Thus, our results on Lao.5Sri.5Mn04 
are consistent with the suggestion in the context of 
LaMnOa that the JT distortion to be not only the driving 
force for the orbital ordering, but also responsible for 
the specific experimental effect of the pronounced tensor 
character of the form factor via the splitting of the 4p 
partial DOS. 

Strong enhancement in the intensity at the energy cor- 
responding to the absorption edge of Mn'^^ was observed 
at the CO superlattice reflection. [Q These results sug- 
gested the presence of a second Mn atom in the system, 
with its absorption edge 4 eV higher than that for the 
Mn'^+ atoms. As the absorption edge for Mn^+ atoms 
lies ~4 eV above the absorption edge for Mn'^+ atoms, 
they claimed that this was a direct evidence of an order- 
ing of the two charge species - Mn'^+ and Mn'*+ in this 
system. From the present calculations, we find that as 
the environment for the Mn^+-like sites is different from 
the environment for the Mn''+-like sites, there is a sub- 
stantial modification in the Mn Ap PDOS, though the 
net charge associated with these two sites, 4.53 at Mn 
1 and 3, and 4.51 at Mn 2 and 4, are almost identical. 
Hence, the modifications at the so-called CO superlattice 
reflection cannot be interpreted on the basis of a substan- 
tial difference in the charge state of the two inequivalent 
sites. Interestingly, our results for the electronic struc- 
ture at these two sites provide a natural explanation for 
the observed results. With the x- and y-axis defined in 
Fig. 1 we have the 4,px states at Mn'^+(1) at lower energies 
compared to the Mn Ap^ states at Mn^+ site (Fig. 3). As 
a result we have an intensity maximum at the CO super- 
lattice reflection. If the axes are rotated by 45 degrees, 
one finds that the Mn Ap^' states at Mn^+(1) and those 
at the Mn^+ site almost coincide (inset in Fig. 4). This 
results in a minimum in the intensity. A similar angu- 
lar dependence is also observed at the 00 superlattice 
reflection. Our reinterpretation is supported by recent 
experiments on Ndo.sSrp.sMnOs which also found a 
similar angular dependence of the intensity at a CO su- 
perlattice reflection. 

In conclusion, we have carried out an analysis of the ex- 
perimental observation for orbital and charge orderings in 
Lao.5Sri.5Mn04. Our results indicate the presence of two 
Mn species with very different environments. One of the 
Mn species has a JT distortion of the oxygens surround- 



ing that atom and hence, in this sense, behaves like an 
Mn'^^ species. The distortions around the so-called Mn^'^ 
atoms are different from the simple breathing mode dis- 
tortion suggested earlier. Extrapolating the magnitude 
of the JT distortion around the Mn^+ to equal the value 
observed in LaMnOa is necessary to obtain the correct 
description of the ground magnetic state. Further, the 
energy dependence at the orbital ordering as well as the 
charge ordering superlattice reflections can be explained 
by the distortions of the MnOg octahedra. The present 
results establish that no real-space substantial charge or- 
dering needs to be invoked in this system in order to ex- 
plain the experimentally observed results; instead band 
structure effects are responsible for the dependence of 
the superlattice intensities on the energy and the orien- 
tation. The above interpretation is consistent with the 
observation that both charge ordering as well as the or- 
bital ordering show the same temperature dependence. 
As the origin of the intensity variation observed at both 
reflections is the same, namely the lattice distortions, our 
results explains why the charge ordering as well as the or- 
bital ordering transition occur over the same temperature 
range in this system. 
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FIG. 1. (a). A schematic diagram of the spin, charge 
and orbital ordering of the Mn atoms in the ab plane. The 
magnetic structure consists of the zigzag chain (thick solid 
line connecting atoms 1, 2, 3 and 4), coupled antiferromag- 
netically to the neighboring chains, (b). The direction of 
displacements of the oxygen atoms when the atom positions 
are relaxed are shown by the thick arrows, while the displace- 
ments of the Mn*+ atoms are indicated by the small arrows. 
The two choices of coordinate axes x and y and x' and y' used 
are shown. 



FIG. 2. The Mn Ap^ (sohd line) and ipy (dotted hne) 
partial DOS projected onto the Mn^^ site showing £^3^2.^2 
ordering. 



FIG. 3. The Ap^o partial DOS at the Mn^+ (sohd line) 
and Mn*+ sites (dotted line) are shown for the x and y axes 
defined in Fig. 1, and in the inset for x' and y' axes defined 
in Fig. 1. 
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